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We introduce a new paradigm for one-dimensional uniform electron gases (UEGs). In this model, n electrons 
are confined to a ring and interact via a bare Coulomb operator. We use Rayleigh-Schrodinger perturbation 
theory to show that, in the high-density regime, the ground-state reduced (i.e. per electron) energy can be 
expanded as e(rs,n) = €t){n)rj'^ + ei(n)rj^ -I- €2(71.) -I- e^{n)rs -f . . ., where is the Seitz radius. We use 
strong-coupling perturbation theory and show that, in the low-density regime, the reduced energy can be 
expanded as e(rs, n) — 'riQ{n)r^^ + rji[n)rs '^^^ + Tl2{n)r'^'^ + . . .. We report explicit expressions for tQ{n), ei(n), 
^2{n)^ e3(n), TjQ[n) and r]i(n) and derive the thermodynamic (large-n) limits of each of these. Finally, we 
perform numerical studies of UEGs with n — 2,3, ... , 10, using Hylleraas-type and quantum Monte Carlo 
methods, and combine these with the perturbative results to obtain a picture of the behavior of the new model 
over the full range of n and values. 
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I. INTRODUCTION 



i—^ In a recent papei^we showed that the traditional concept of the uniform electron gas (UEG), i.e. a homogeneous 

C/3 system of finite density, consisting of an infinite number of electrons in an infinite volumejl^"^ is inadequate to model 
^ the UEGs that arise in finite systems. Accordingly, we proposed to embark on a comprehensive study of quasi-exact 

properties of finite-size UEGs, in order eventually to create improved approximations in density-functional theory.^ 
^ In an earlier paper, ^ we introduced an alternative paradigm, in which n electrons are confined to a Z)-sphere (with 

i( D > 2), that is, the surface oi a {D + l)-dimensional ball. These systems possess uniform densities, even for finite n 
and, because all points on a ZJ-sphere are equivalent, their mathematical analysis is relatively straightforward.'^'ti^ In 
Q the present paper, we study the one-dimensional (D — 1) version of model, in which n electrons are confined to a ring 
Q of radius R. The electron density of this n-electron UEG, which wc will call n-ringium, is 

n I , 

where = nR/n is the Seitz radius. In this study, the high-density (small-rs) limit is defined by i? — )■ for fixed n, 
\^ while the low-density (large- r^) limit is defined by i? — )■ (X) for fixed n. We do not include a fictitious uniform positive 
\^ background charge because, unlike the situation in 2D and 3D UEGs, its inclusion in ID systems causes the Coulomb 

energy to diverge. 

In most previous work on the one-dimensional (ID) UEG, the true Coulomb potential l/ri2 has been avoided 
because of the intractability of its Fourier transform. Instead, most workers have s oftened t he potential, either by 
T— I adding a transverse harmonic componenlli^HIIl or by using a potential of the form l/yVf^~T7t^. In the latter case, the 
/i parameter eliminates the singularity at ri2 = while retaining the long-range Coulomb tail.E^lti^ 

However, the introduction of a parameter /i > is undesirable, for it modifies the physics of the system in the 
high-density regime where neighboring electrons repel far too weakly. It is also unnecessary, because the true Coulomb 
^ potential is so repulsive that it causes the wave function to vanish when any two electrons touch, thereby removing the 
possibility of an energy divergence.^ For ID systems, we have recently shown that the exact wave function ^' behaves 
as 

*(^i2)=ri2(l + ^)+0(rL) (2) 
for small ri2,'^ which is the ID analog of the (three-dimensional) Kato cusp condition .1221 
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This nodal behavior leads to the ID Bose-Fermi mappin^^ which states that the ground state wave function of the 
bosonic (B) and ferniionic (F) states are related by 5'b(R) = |^f(R-)|, where R = (ri,r2, . . . ,r„) are the one-particle 
coordinates. In case of bosons, the divergence of the Coulomb potential has the effect of mimicking the Pauli principle 
which prohibits two fermions from touching. This implies that, for ID systems, the bosonic and fermionic ground 
states are degenerate and the system is "spin-blind" . Consequently, the paramagnetic and ferromagnetic states are 
degenerate and we will consider only the latter.^i^ 

The electrons-on-a-ring paradigm has been intensively studied as a model for quantum rings (QRs), which are tiny, 
self-organised, ring-shaped semiconductors^^ ■^'^ characterised by three parameters: radius (R), width (S) and electron 
number (n). M odern microfabrication technology has yielded InGaAs and GaAlAs/GaAs QRs that bind only a few 
electronsp'^Il jn contrast with the mesoscopic rings on GaAs which hold much larger numbers of electrons .'^S' These 
low-dimensional systems are the su bject of considerable scientific interest and have been intensively studie d, bo th 
experimentall5^22H23 and theoretically,'i^E2lMIlIlniainly because of the observation of Aharonov-Bohm oscillations .'^^'^^Hlll 

As a first approximation, QRs can be modelled by electrons confined to a perfect ring (i.e. 6 = 0). In a recent 
paper^ we considered a pair of electrons (i.e. 71 = 2) on such a ring and discovered that their Schrodinger equation 
can be solved exactly, provided that the radius takes one of an infinite number of special values. Some of the solutions 
exhibit the Berry phase phenomenon, i.e. if one of the electrons moves once around the ring and returns to its starting 
point, the wave function of the system changes sign. QRs are among the simplest systems with this peculiar property. 

In Section [nj we first use Ray leigh- Schrodinger perturbation theory to investigate the energy in the high-density 
regime^'' and then strong-co upli ng perturbation theory to study the low-density regime, where the electrons form a 
Wigner crystal.l^In Section [llH we use explicitly correlated (EC) methods to determine the energy of n-ringium for 
n = 2, 3, 4, 5. These methods are accurate for very small n but their cost grows very rapidly with n. In Section IV we 
turn to quantum Monte Carlo (QMC) approaches for studying n-ringium up to n = 10 over a range of densities. These 
methods provide a different approach to the many-body problem: variational Monte Carlo (YMCj^^^^ and diffusion 
Monte Carlo (DMC)AS..Sfli methods can be used to treat systems in one and higher dimensions at a computational cost 
that grows relatively slowly with n (at least, when n is not too large.'^ 

We frame our discussion in terms of reduced energy e(rs,n), i.e. energy per electron, so that we can pass smoothly 
from finite to infinite n. One of the key goals of the paper is to develop an understanding of the correlation energy, 
which is defined as the difference 

ecirs,n) = e{rs,n) ~ e-iiFirs,n), (3) 

between the exact and Hartree-Fock (HF) energies. Atomic units are used throughout, but we report total energies in 
hartrees (-Eh) a-nd correlation energies in millihartrees (mE\^). 

II. PERTURBATIVE METHODS 
A. High-density expansion 

The Hamiltonian of the system is 

1 " 92 " 1 



i — 1 *■ i<j 

where 9i is the angle of electron i around the ring center, and 



= |r, - I = R^2 - 2 cos(0, - 9^) (5) 

is the across-the-ring distance between electrons i and j. 

In the high-density [i.e. small Ts) regime, the kinetic energy is dominant and it is natural to define a zeroth-order 
Hamiltonian 

and a perturbation 
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The non-interacting orbitals and orbital energies are 

XaW - (2^i?)"'/'exp(m0), (8) 

«a = ^, (9) 

where 

[...,-2,-1,0, +1, +2,..., if n is odd, 
a — < ■? 1 1 ^ 

I . . . , -|, -|, +|, +|, . . . , if n is even. 

A Slater determinant "4*2 of any n of these orbitals has an energy Ei and is an antisymmetric eigenfunction of Hq. In 
the lowest energy (aufbau) determinant "4'o, we occupy the orbitals with 

n— Iri — 3 n — 3ri — 1 , ^ 

a = , ,H ■ (11) 

2 2 2 2 ^ ' 

Following the approach of Mitas,'^ one discovers the remarkable result 

n 

"*o«n^.j-, (12) 

where 

'e,-e 



2flsin \ M (13) 



is a signed interelectronic distance. It follows immediately that "V&o has a node whenever 0i = 9j and, therefore, 
possesses the same nodes as the exact wave function. This will have important ramifications in Section |IV[ 
Rayleigh-Schrodinger theory yields the perturbation expansion for the reduced energy 

e(r„ n) = + + €2{n) + e^{n)r, + . . . , (14) 

where the high-density coefhcients (n) are found by setting R — 1 and evaluating 

^2 



eo{n) = ^(l-olFol^o), (15a) 



ii{n) = -!L{4.„|V-|*„), (15b) 



(150 



Eq " Ei 



^ {Eq - Ei){EQ — Ej) 



1. Double-bar integrals 

To evaluate the coefficients ej{n) with j > 0, one requires the "double-bar" integrals 

iabWcd) . rr X:(^i)x^(^.)[Xe(^0x.(^2)-X.(^2)x.(^0],,^,,^^ 
Jo Jo ri2 
By elementary integration, one can show that 

(afc||ccJ) = < „ . (17) 

0, otherwise. 



where 



Vp,,= ^[i;{p+^)-^P{q+l)], (18) 



and is the digamma function.lS^ 
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TABLE I. High-density coefficients for n-ringium. (("(3) is Apery's constant!^ 



n 


eo(n) 


ei(n) 


e2(n) 




2 


1 ^2 
32 


1 

2 


1- ^ 


8(121n2-19) 1 16{26-7C(3)) 


3 


1 ^2 
27^ 


20 
27 


16 1436 
9 8l7r2 


8(1080 ln2-997) , 8(13046-4725C(3)) 
81?r2 ' 2437r'> 


4 


5 2 
128 


9 

10 


109 244168 
45 10125ir^ 


0.00487354 


5 


1 ^2 
25^ 


892 
875 


4688 514012364 
1575 173643757r2 


0.00556461 


6 


_35_ 2 
864 " 


6323 
5670 


2339 461265158 
675 133953757r^ 


0.00605813 


7 


2 ^2 
49 


13528 
11319 


1420256 33870168846728 
363825 8736329621257r^ 


0.00642454 


8 


21 ^2 
512 


7591 
6006 


20349053 81975019672689056 
4729725 19193716177886257r2 


0.00670533 


Q 


10 _.2 

243 


4831544 
3648645 


66244064 266761139809046216 
14189175 57581148533658757r^ 




10 


33 2 
800 " 


2512297 


1207979879 7026989855398034506022 


0.00710359 


1823250 


241215975 1414480913729327193757r2 


oo 


7rV24 


In ^/n 


-7rV360 


0.00844621 


2. 


Zeroth order 









The zeroth-order coefficient (15a) becomes 



9 OCC 9 

TT"^ \ - a 



(19) 



where the "occ" indicates sums over all occupied orbitals (11), and this reduces to 

, , n2 - 1 TT^ 



n2 24 



In the thermodynamic (i.e. n ^ oo) limit, this approaches 



TT 

24' 



which is identical to the kinetic energy coefficient in the ideal Fermi gas in ID. fulfil 



3. First order 



(20) 



(21) 



The first-order coefficient (15b) becomes 



which can be reduced to 



ei{n) 



eiin) = — V(a6||a6), 

7?/^ ^ ^ 



a<b 



(22) 



(23) 



This can be found in closed form for any n (see Table |T|. Because of the slow decay of the Coulomb operator, the 
coefficient grows logarithmically with n and it can be shown that 



ei(n) - In ^^:+ (In 2 + 7/2 -3/4) + 0(^-2 Inn), 

where 7 is the Euler-Mascheroni constant. 

The sum of the first two terms in ( 14 ) gives the HF energy of n-ringium 



/ N <^o{n) ei n) 
eiiF[rs,n) = — — H . 



(24) 



(25) 
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Second order 



The second-order coefficient ( 15c I becomes 

= (a6||r.)(r.||a6) ^^6) 

^ a<b r<s '^^ ~^ '^'^ 

where the "virt" indicates sums over ah virtual orbitals. If the double-bar integrals do not vanish, i.e. a -f 6 = r -f s, 
then 

Kr + Ks- Ka- {r - a){r -b), (27) 

and we obtain 

_j OCC C30 T^2 

e2(n) = -iV V . (28) 
^ ^ (r-a)(r-b) 

a<b r=r,„i„ ^ ' ^ ' 

where 

n -\- 1 

r^in = \- max(a + 6,0). (29) 

The sums in (28 1 can be evaluated in closed form for any n (see Table [l]). In the — ?> limit, the higher terms in (14) 
vanish and the e2{n) expressions in Table |l] are therefore the exact correlation energies of infinitely dense n-ringium. 
In the thermodynamic limit, e2("-) approaches 

Uln 

7T \ r—b 



OCC oo 



£2 = - hm - 



n->oo n ^ ^ (r — a)(r — b) 

a<br=r„-,i„ ^ ' 

— ) dxdy 



which implies that, in the dual thermodynamic/high-density limit, the exact correlation energy of ringium is —27.4 mE\^ 
per electron. The same value of £2 can be derived for ID jellium,'^ affirming the equivalence of the electrons-on-a-ring 
and electrons-on-a-wire models in the thermodynamic limit.'^ 

Using a quasi-lD model with a transverse harmonic potential, Casula et al. were led to conclude that, in the same 
limit, the correlation energy vanishes.^ This qualitatively different prediction stresses the importance of employing a 
realistic Coulomb operator for high-density UEGs. 



5. Third order 



The third-order coefficient (15d) becomes 



1 {ab\\r s) {r s\\cd) {cd\\ab) 1 {ab\\rs){rs\\tu){tu\\ab) 

^ Z-^ (r - n\(r - h\(r - A(t - rl\ + 8^ 2^ 



1 {ab\\rs){cs\\tb){rt\\ac) 1 {ab\\rs){ar\\ct){rs\\ab) , , 



abc rst 



{r - a){r - b){r - a){r - c) ti" ^ {r - a){r ~ b){r - a){r - b) 



and, like £2(^1), this can be rewritten in terms of products of Vp q. The expression is cumbersome but can be evaluated 
in closed form for any n and Table |T] illustrates this for n = 2 and 3. 
In the thermodynamic limit, £3(n) approaches the numerical value 

£3 = +0.00844621, (32) 
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but we have been unable to obtain this in closed form. Numerical evidence suggestgSf that (32) is also true of ID 
jellium. 

Interestingly, second- and third-order perturbation theories applied to ID jellium do not encou nter d ivergence issues 
as in 2D and 3D jellium, where one has to use resummation techniques to produce finite results !^2155j ^j^g divergence 
occurs from third order and second order for 2D jellium and 3D jellium, respectively. In the case of ID jellium, every 
terms of the perturbation expansion seem to converge. 



B. Low-density expansion 



In the low-density (r^ > 2) regime,'^ the electrons form a Wigner crystal. Using strong-coupling perturbation 
theory,'^ the energy can be written 



3/2 



(33) 



where the first term represents the classical Coulomb energy of the static electrons and the second is their harmonic 
zero-point vibrational energy. 

The Wigner crystal, which is the solution to the ID Thomson problem,!^ consists of n electrons separated by an 
angle 27r/n and yields 



n-l 



2n2 sin(fc7r/n) 



(34) 



The second term in the expansion (33) is found by summing the frequencies of the normal modes obtained by 
diagonalization of the Hessian matrix. For electrons on a ring, the Hessian is circulant and its eigenvalues and 
eigenvectors can be found in compact form, yielding 



771 (n) 



r3/2 



2 — sin^(A:7r/n) . 



4„5/2 E ^ E ^ sin3(VJ)n)"' 

1=1 \ k=l ^ ' ' 



(35) 



In the thermodynamic limit, one finds that 



, ^ ln(2/7r) +7 , 
7?o = In x/r^ + ' ^ '- + o(n°), 



(36) 



which has the same logarithmic divergence as ei, but with a different constant term. Likewise, one can show that 



r3/2 



E 



4^5/2 ^ ^ {kn/nf 
^ ^2Li3(l)-Li3(e^«)-Li3(e-*«)d0, 



■ sin^(ifc7r/n) 



(37) 



where Lis is the trilogarithm function.!^ We have not been able to find this integral in closed form, but it can be 
computed numerically with high precision, and yields ryi = 0. 359933, w hich is identical to the value found by Fogleil^ 
for an infinite ultrathin wire and a potential of the form 1/ + M^- This shows that, unlike the high-density limit 
where the details of the interelectronic potential are critically important, the correct low-density result can be obtained 
by using a softened Coulomb potential. 

Thus, in the dual thermodynamic/low-density region, we have 



ln(V27r) - 3/4 0.359933 



0(r, 



(38) 



The same expansion can be derived for the infinite wire^ confirming the equivalence of the electrons-on-a-ring and 
electrons-on-a-wire models in the thermodynamic limit .121 
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TABLE II. Convergence with M of the energy of 2-ringium with rs — 1. 



M £(1,2) -ec(l,2) 

0.808 425137 534 

1 0.797 201143 955 11.223 993 579 

2 0.797175 502 306 11.249 635 229 

3 0.797175 223 852 11.249 913 682 

4 0.797175 219 345 11.249 918190 

5 0.797175 219 257 11.249 918 277 

6 0.797175 219 255 11.249 918 279 



III. EXPLICITLY CORRELATED METHODS 



Because the full set of interelectronic distances determine the positions of the electrons to within an overall 
rotation that is irrelevant in the ground state, it is appropriate to adopt these variables as natural coordinates and to 
expand the correlated wave function in terms of these distances. 



A. 2-ringium 



The HF wave function for 2-ringium is 

^^o = fi2- (39) 

In the light of its simplicity, and following our previous analysis of the quasi-exact solutions,"^ it is natural to consider 
correlated wave functions that are products of ^'I'o and a correlation factor, viz. 

M 

2 



*A/ = ^'"''"s- (40) 



The overlap, kinetic and potential matrix elements can be found as outlined in the Appendix. 

Table [n] shows the energies obtained by solving the secular eigenvalue problem for Ts = 1. They converge rapidly, 
with M = 1, 2, 4, 6 yielding milli-, micro-, nano- and pico-hartree accuracy, respectively. It is interesting to compare 
the correlation energy (—11 mii^h) with the corresponding value (—114 mi?h) for two electrons on a 2D sphere.'^'' One 
normally expects the correlation energy to decrease in higher dimensions^® but the ID case is anomalous because the 



HF wavefunction ( 12 ) places the two electrons in different orbitals. 



The key discovery from this investigation is that including just the linear (ri2) and quadratic (r^2) terms in the 



expansion ( 40 ) affords microhartree accuracy for the energy of 2-ringium. We now ask whether this is true for larger 
values of n. 



B. 3-ringium 

The HF wave function for 3-ringium is 

^^-o = fi2?'l3 7'23, (41) 

and we have explored both Hylleraas-type wavefunction^S^HUl 

hi 

'*m' = '*oE E c.,fesl4s3^ (42a) 

7Ti=0 i+2j+3fe<m 

Si = r-12 + r-ia + r23, (42b) 

S2 ?'12»'13 +?'12T'23 +fl3?'23, (42c) 
S3 = fl2''l3?'23, (42d) 
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TABLE III. Convergence with M of the energy of 3-ringium with Vs — 1. 



M 


HyUeraas expansion 


Jastrow 


expansion 


e^i, 6) 


— Eel,!, O) 


6(1,3) 


-ec(l,3) 





1.106 281644 485 





1.106 281644485 





1 


1.091649 204 702 


14.632 439 783 


1.090 999 267 912 


15.282 376 573 


2 


1.090 936176 037 


15.345 468 448 


1.090 936 808 374 


15.344 836111 


3 


1.090 935619110 


15.346 025 375 


1.090 936 772 712 


15.344 871 773 


4 


1.090 935 608 007 


15.346 036 478 


1.090 936 607 858 


15.345 036 627 


5 


1.090 935 607817 


15.346 036 667 


1.090 936 593 657 


15.345 050 828 


6 


1.090 935 607811 


15.346 036 674 


1.090 936 589183 


15.345 055 301 


7 


1.090 935 607810 


15.346 036 674 


1.090 936 588 261 


15.345 056 224 





TABLE IV. 


Convergence with Af of the energies of 4- 


and 5-ringium with rs — I. 




M 




4-ringium 


5-ringium 




^(1,4) 


-fc(l,4) 


e(l,5) 


-ec(l,5) 





1.285 531 





1.414213 





1 


1.269 785 


15.746 


1.398192 


16.021 


2 


1.268 259 


17.272 







and Jastrow-type wavefunctioni 



1621 



M 



E 





(43) 



The required matrix elements can be found as outlined in the Appendix. 

The HyUeraas expansion converges rapidly for 3-ringium with = 1 and Table III reveals that, as in 2-ringium, M 
= 1, 2, 4, 6 yields milli-, micro-, nano- and pico-hartree accuracies, respectively. The reduced correlation energy is 
roughly 35% greater than that in 2-ringium. Because of its factorized form, the limiting Jastrow energy is « 1 /i-Eh 
above the exact value. 



C. 4- and 5-ringium 



The HF wave functions for 4- and 5-ringium, respectively, are 

''^'O = fi2 fi3 fi4 7=23 f24 ^4, 
= ri2 fi3 fi4 ri5 f23 ^24 ?^25 ^34 hb ^^45- 



(44) 
(45) 



HyUeraas calculations on these systems are complicated because of the large number of many-electron integrals which 
are required. Nonetheless, we were able to perform such calculations, up to Af = 2 for 4-ringium and up to M = 1 for 
5-ringium, and the results are summarized in Table IV It is important to allow rijVki terms (which couple two electron 



pairs) and VijVik terms (which describe three-electron interactions) to have distinct HyUeraas coefficients: failing to do 
so raises the energy by sa 1 /ii?h- The energies in Table IV are higher than our best estimates (see Table [v]) by roughly 
2 uiEh (for M = 1) and 50 fiEh (for M = 2). 



IV. QUANTUM MONTE CARLO METHODS 



A. Variational Monte Carlo 



In the VMC method, the expectation value of the Hamiltonian with respect to a trial wave function is obtained 
using a stochastic integration technique. Within this approach a variational trial wave function 4't(R, c) is introduced, 
where c — (ci, C2, . . . , cm) are variational parameters. One then minimizes the energy 

1 /vI'T(R,c)gM/T(R,c)dR 

'^^^ - n /*T(R,c)2dR ' 
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with respect to the parameters c using the Metropolis Monte Carlo method of integration.^^ The resulting VMC 
energy is an upper bound to the exact ground-state energy, within the Monte Carlo error. Unfortunately, any resulting 
observables are biased by the form of the trial wave function, and the method is therefore only as good as the chosen 
\['t. 

Here, we use electron-by-electron sampling with a transition probability density given by a Gaussian centered on the 
initial electron position. The VMC time step, which is the variance of the transition probability, is chosen to achieve a 
50% acceptance ratio.l^ 



B. DifFusion Monte Carlo 

DMC is a stochastic projector technique for solving the many-body Schrodinger equation.l^^Hini its starting point is 
the time-dependent Schrodinger equation in imaginary time 

' =iH-S)^{R,r), (47) 



and it is exact, within statistical errors. For r — >■ oo, the steady-state solution of Eq. (47 1 for S close to the ground-state 
energy is the ground-state \1/(R).'^DMC generates configurations distributed according to the product of the trial and 
exact ground-state wave functions. If the trial wave function has the correct nodes, the DMC method yields the exact 
energy, within a statistical error that can be made arbitrarily small by increasing the number of Mon te Carlo steps. 
Thus, as in VMC, a high quality trial wave function is essential in order to achieve high accuracy.'^^'^ 

Our DMC code follows the implementation of Reynolds et al.^ using a population of ^ 5000 walkers for each 
calculation. We have carefully checked that the population-control bias is negligible. The dependence of the energy 
upon the DMC time step At was also investigated and the extrapolated value of the energy at Ar = is obtained by 
a linear extrapolation. The number of points used in the fitting procedure depends on . A minimum of 4 points has 
been used for linear interpolation in the set Ar = 0.0001, 0.0002, 0.0005, 0.001, 0.002 and 0.005. The extrapolated 
standard error is obtained by assuming that the data follow a Gaussian distribution.^ We note that the algorithm 
developed in Ref. [65! does not significantly reduce the time-step error in the present case. 



C. Trial wave functions 



We have employed Jastrow trial wave functions 

n M 

"*M = "*onE^fc^S' (48) 

i<j k—1 

choosing M — 5 in order to obtain microhartree energy accuracy for rs — 1. The co efficien ts Ck were optimized using 
Newton's method following the methodology developed by Umrigar and co- workers. 1^3211 Por < 1, we used energy 
minimization; for > 1, energy minimization was unstable and we minimized the variance of the local energy.!^ 



D. Fixed-node approximation 



DMC algorithms can be frustrated by the sign problem in fermionic systems.'S^l^ To avoid this, it is common to 
apply the fixed-node approximation, i.e to write the wave function as the product of a non- negative function and 
a function with a fixed nodal surface. The DMC method then finds the best energy for that chosen nodal surface, 
providing an upper bound for the ground-state energy. The exact ground-state energy is reached only if the nodal 
surface is exact but, fortunately for us, the nodal surface of the HF wave function Eq. (12 1 is exact and, therefore, 
DMC calculations using the trial wave function ( 48 1 yield the exact energy. We have no fixed- node error. 



E. Results and discussion 



Table [V| summarizes the results of a systematic study of n-ringium systems with rg = 1. In all cases, our DMC 
calculations yielded energies with statistical uncertainties within 1 /ii?h and this allowed us to assess the accuracies of 
our explicitly correlated calculations. 
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TABLE V. Hartree-Fock, Explicitly Correlated and Diffusion Monte Carlo energies of n-ringium with Vs — 1. Statistical errors 
in the last digit of the DMC energies are shown in parentheses. 



n 


eHF(l, n) 


eEc(l,n) 


eDMc(l,'"-) 


2 


0.808 425 


0.797175 


0.797175(0) 


o 


1 infi9>?9 

-L . -LUU ZOi 


1 .uyu you 




4 


1.285 531 


1.268 259 


1.268212(1) 


5 


1.414213 


1.398192 


1.395 774(1) 


6 


1.514 978 




1.495 841(1) 


7 


1.598 000 




1.578393(1) 


8 


1.668 711 




1.648 770(1) 


9 


1.730 359 




1.710172(1) 


10 


1.785 044 




1.764 671(1) 




5 10 15 20 

r, (a.u.) 



FIG. 1. Reduced correlation energies (in m_Eh) for n-ringium with various Vs- 



Table VI summarizes our best estimates of the correlation energies of n-ringium for various (see also Fig. [T]) . For 



— 0, we use the exact e2{n) values from Table |lj For rg — 0.1, we use the Fade approximant 

ecirs,n)^ , (49) 

1 - [e3{n)/e2{n)\rs 

which provides microhartree accuracy. For n = 2 and n = 3, we use the Expli citly Correlated results from Table \V\ 
For n = oo and 1 < < 20, we use the DMC results from Lee and Drummond.'i^ For n = oo and 0.2 < < 0.5, we 
performed DMC calculations using the CASINO softward^ following the Lee-Drummond methodology. 

For the remaining cases (4 < n < 10 and Vg > 0.2), we used our own DMC program. We achieve su.h-fiE\^ 
uncertainties for > 1 (where the electrons become localized and approach a Wigner crystal^ but it is difficult to 
achieve this for smaller r^, where the uncertainties are 10 - 40 /J.i?h- 



V. CONCLUSIONS 



We have studied n-ringium using explicitly correlated and quantum Monte Carlo methods. Using Hylleraas wave 
functions, we have obtained the near-exact ground-state energy of the n = 2 and n = 3 systems for various values 
of the Seitz radius Ts ■ For n > 4, we have performed exact-node DMC calculations to find the exact ground-state 
energies, with statistical errors in the fiEj^ range. 

We have shown that the reduced correlation energy of n-ringium is 

ec{rs,n) ^ €2(11) + e3{n)rs + . . . (50) 

for high densities, and 



(51) 
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TABLE VI. Reduced correlation energies (in miJh) for n-ringium with various r^. Statistical errors in the last digit of the DMC 
energies are shown in parentheses. 



Ti\Ts 





0.1 


0.2 


0.5 


1 


5 


10 


20 


2 


13.212 


12.985 


12.766 


12.152 


11.250 


7.111 


4.938 


3.122 


3 


18.484 


18.107 


17.747 


16.755 


15.346 


9.369 


6.427 


4.029 


4 


21.174 


20.698 


20.24(2) 


19.00(1) 


17.320(1) 


10.390(0) 


7.085(0) 


4.425(0) 


5 


22.756 


22.213 


21.66(2) 


20.33(1) 


18.439(1) 


10.946(0) 


7.439(0) 


4.636(0) 


6 


23.775 


23.184 


22.63(2) 


21.14(1) 


19.137(1) 


11.285(0) 


7.653(0) 


4.762(0) 


7 


24.476 


23.850 


23.24(2) 


21.70(1) 


19.607(1) 


11.509(0) 


7.795(0) 


4.844(0) 


8 


24.981 


24.328 


23.69(3) 


22.11(1) 


19.940(1) 


11.664(0) 


7.890(0) 


4.901(0) 


9 


25.360 


24.686 


24.04(2) 


22.39(1) 


20.186(1) 


11.777(0) 


7.960(0) 


4.941(0) 


10 


25.651 


24.960 


24.25(4) 


22.62(1) 


20.373(1) 


11.857(0) 


8.013(0) 


4.973(0) 



27.416 26.597 25.91(1) 23.962(1) 21.444(0) 12.318(0) 8.292(0) 5.133(0) 



for low densities. Expressions for the coefficients are given in Eqs. (28 1, (31 1, (34), (24) and (351 
In the thermodynamic limit, we have found that 



— +0.008446 r, + ..., 
360 

h-i{V2n) - 3/4 0.359933 

^ ^ 3/2 



(52) 
(53) 



and shown that the ringium and jellium models are equivalent in the thermodynamic limit. 

This provides a detailed picture of the energy of this new model over a wide range of n and values and we believe 
that the correlation energies in Table VI are the most accurate yet reported for rt-ringium. These systems are distinct 
uniform electron gased^ and can be used to design a new correlation functional for ID systems. We will report such a 
fimctional in a forthcoming paper.^^ 
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APPENDIX 



For the ground state, the Hamiltonian ^ can be recast as 



i<3 



4i?2 



drl ' AR^dr. 



r,j d 



n n 2 , 2 _ 

' ■•' ^ ' jk ' ij ' ij 



EE- 



52 



2 ^ik'^jk 



(54) 



The first term in (54 1 contains the two-body parts of the Hamiltonian while the second includes coupling between 
electron pairs. 

The n-electron overlap integrals needed in calculations on 7i-ringium can be systematically constructed using the 
unit-ring Fourier resolution 



where 



Br, 



(-l)'=m! 



(m/2 + A:)!(TO/2-fc)! 



(55) 



(56) 
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is a signed binomial coefficient. Eq. ( 55 ) is valid for m > and terminates if m is an even integer. 

Resolving each integrand factor, swapping the order of summation and integration, performing the integrations and 
resumming, often leads to beautiful expressions. For example, the cyclic n-electron integral yields 



a b c 
''12^23^34 ■ 



''n,l) — ^ Ba^kBb,kBcM ■ ■ ■ Bz,k, (57) 



fc— — oo 



which can be written as a n+iFn hypergeometric function of unit argument.'^ 
In some cases, the sums can be found in closed form, for example. 



= Sa.o, (58) 
{r'^Ai) - Ba.oBbfl, (59) 



a\h\c\ ( 


a+b+c ^ 
2 ) 


! 


(^a+b\ 


! 


( a+c\ 
\ 2 ) 




!(t; 




!(§) 


! 



(^?2^23^3l) = / TTT /.^ N / ^ /\ ■ (60) 



but this is not possible in generahf^ (Eq. (60) uses the result of Problem 62 of Knuth's book.'^ However, sums such 



as (57) converge rapidly and are numerically satisfactory. 
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